Particle Trajectories in Motel Current Sheets, 
Part I: Analytical Solutions 


N66-18363 


(ACCESSION NUMBER) 




i CR OR TMX OR AD NUMBER) 


n 


by 

T. W. Speiser* 


(CATEGORY) 




GPO PRICE $. 
CFSTI PRICE(S) S . 

Hard copy (HC) 
Microfiche (MF) 

ff 653 July 65 


Abstract . Approximate analytical solutions are found for two 
model current sheets . In the first the magnetic field is linear and 
reverses across a neutral sheet, and the electric field is everywhere 
uniform, perpendicular to the magnetic field and parallel to the neutral 
sheet. Charged particles of either sign never come out of the neutral 
sheet and their energies increase without bound. In the second model a 
small component of the magnetic field perpendicular to the neutral sheet 
is added. This component not only serves to bring particles out of the 
sheet, but turns protons and electrons toward the same direction, 9 C 0 away 
from the accelerating electric field. The particles are accelerated and 
then ejected when they have been turned $ 0 °, and. the emergent pitch angles 
to a magnetic line of force will be small if the perpendicular magnetic 
field component is small. 


Introduction . Hess [1964] has reported the discovery of a magneticall 
neutral sheet in the earth' s geomagnetic tail, with measurements taken on 
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board the IMP-1 satellite. A somewhat uniform magnetic field of about 
20 y ( 1 y — 10 ^ gauss = 10 ^ weber/m 2 ) is found from about 10 R (earth- 
radii) to at least 30 R g in the anti-solar direction. This field ^.s 
predominantly in the solar direction above a plane (roughly identified 
as the magnetic equatorial plane) , and in the anti-solar direction below 
this plane. The magnetic field reverses across a sheet of thickness 
about 0.1 R g and goes to zero within this neutral sheet. 

Ocher neutral-, or more generally current-, sheets may occur in 
other situations, such as a day- side magnetospheric current sheet, 
neutral sheets occuring in interplanetary space, neutral sheets associated 
with solar flares, etc. 

It is of interest to look at charged particle trajectories about 
such sheets. Adiabatic theory cannot be used across such a neutral sheet, 
because the magnetic field changes significantly in distances much less 
than a gyroradius. The charged particle equations of motion must there- 
fore be either solved analytically, or computed numerically. 

Thi s paper is concerned with analytical solutions ir, two model 
current sheets. Part II { Spei ser , 19 o 5~] applies these results to a 
magnetospheric tail model and also discusses numerical results. 

A simple linear model 

The simplest model for the fields about a neutral sheet which can 
be discussed analytically is: 

B = - b _ e , 
d V 
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a e 


( 2 ) 


E = - 
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where b is the strength of the magnetic field when x = d, the sheet half- 
thickness, and a is the strength of the electric field. The physical 
significance of such an electric field will be discussed in Part II 
f Spei ser 9651 • This field will merely be assumed for the present 
treatment . 

The co-ordinate system being used is sketched in Figure 1 . The 
above magnetic and electric fields are also indicated in Figure 1 . 

The equations of motion for a particle in the neutral sheet, using 
these fields are 


x = C i z x 

(3) 

y = o 

(4) 

z = - C - C , x x 

(5) 


q b q 

where C, = , and C, = - a. Mks units are used. 

, 1 m d 3 m 

Speiser [1964, (a), (b)] has given the solution to these equations 
for large time. (See Appendix A.) The result is that the particle 
executes a damped oscillation about x = 0 (the amplitude going as 
1 / 1 1 Z' 4 ) , while accelerating in the +e z direction for electrons and the 
-e direction for protons. This can be understood as follows. The 
electric field in equation (5) accelerates a proton in the -e z 


direction . 


z then becomes proportional to -t ignoring the term -C 1 xx 
for the moment. Equation (3) then becomes 

t 

x = - k x 

where k = j zj , or k = ^ ' — t * 

m d 

Equation (6) is just the equation for oscillatory motion, wirh 
spring constant k. k gets larger with time, however, implying that the 
spring gets stiffer with time, thus the amplitude of oscillation decays. 

The term -Cix x in equation (5) is thus oscillatory and approaches 
zero, justifying its neglect but the term is not neglected initially 
in Appendix A. k is proportional to q 2 , so particles of either sign 
will oscillate. 

The oscillation in x(t) is due to the V x B force which is always 
toward x = 0 for x either positive or negative, because of the reversal 
of the magnetic field. 

The net result of this simple model is that charged particles of 
either sign never come out of the neutral sheet, and their energy in- 
creases without bound. See Figure 2 for a sketch of the results of 
this simple model. 

The linear model with small perpendicular field added . Consider 
now the addition of a small magnetic field component perpendicular 
to the neutral sheet, i.e., in the + e x direction referring to the 
co-ordinate system of Figure 1. The magnetic field of equation (l) 


now becomes: 


r 



where T| will be assumed small. 

The equations of motion become: 

x = C i x z ’ ( 3 ) 

y = C z T] i (9) 

z = - C 3 - Ci x x - C 2 Tj y (10) 

where C ^ = (q/m) (b/d) , C 2 = (q/m) b, C 3 = (q/m)a, C, and C 3 being 

the same as for the simple model. When 7] = 0 , the simple model results. 

Even without solving these equations, the particle motion can be 
understood qualitatively as was done for the simple model. 

Consider a proton (the arguments also hold for electrons with 
appropriate changes in sign) incident on this neutral sheet with small 
velocity. (Strictly speaking, the sheet is not now a neutral sheet.) 

The proton will be accelerated initially in the negative z direction 
(see Equation (10)), gaining a velocity, z, proportional to -t as for 
the simple model. From Equation (9) , there will then be an acceleration 
in the -y direction proportional to z or -t, and thus y will be propor- 
tional to -t 2 . (Note that y a C 2 C 3 a (q/m) 2 , so either protons or 
electrons are turned toward the -ey direction.) 

As long as z is negative, Equation (8) will imply oscillatory 
motion in x(t) by the same arguments in the previous section on the 
simple model. Thus, as long as z is negative, the term -C ; x x in 



Equation (lO) is oscillatory and will be assumed small as a first order 
approximation. The oscillations may now be either damped or growing 
depending on whether z increases or decreases with time. The last term 
in Equation (10) grows as + t 2 so that z will grow negatively until z 
goes to zero, and z will then diminish in absolute value going to zero 
and even becoming positive after some time. Thus x(t) will execute 
damped oscillatory motion until z = 0 (the spring gets stiffer) . After 
z becomes positive and until z = 0, x(t) will execute growing oscillatory 
motion. After z goes positive, however, x(t) will no longer oscillate 
but will increase exponentially, thus ejecting the particle from the 
neutral sheet. See Figure 3 for a sketch of the particle trajectories 
in this model. 

Results of the approximate theory . The detailed calculations are 
contained in Appendix B, to which the following alphabetic equation 
indices refer. First integrals of the equations of motion (Equations 
(8), (9) and (10)) are obtained exactly (Equations (i), (j) and (n)). 

A first approximation is used to obtain the time of ejection, t, which 
turns out to be inversely proportional to ?| b (q/m) , (Equation (s)). 
Electrons are therefore ejected from the neutral sheet much sooner 
than protons. The velocity in the -e y direction at the time of ejection 
is found to be independent of (q/m), (Equation (w)), and the maximum 
pitch angle of the emergent particles (Equation (z)) is proportional 
to Tj, and to |xq/u - c| , where u is the bulk flow velocity exterior 
to the neutral sheet, and c is a number which is not very Different 
from 1. This result implies that a ~ 0 for T\ - 0, which means that 



the particles would never come out of the neutral sheet if the small 
perpendicular field were not added. This agrees with the result of the 
simple model, and implies that a will he small if T] is small. Since a is 
also proportional uo jx o/u - cj, particles incident on the current sheet 
with x 0 se u will all emerge with pitch angles close to zero. 

Dur.gey (ljop) suggested that the electric field can be transformed 
away in this model. Calculation using such a transformation are made in 
Appendix C. The particle motion in the moving (transformed} system is 
easily visualized, since only a magnetic field is involved. The perpendicular 
magnetic field component causes a circular drift in the current sheet while 
the x- coordinate is oscillating due to the magnetic field reversal (b ) . 

y 

(See, for example Equations (ac), (ah), and (al) .) The ejection time is 
seen to be just a half-period of the circular drift (Equation (am)) and 
corresponds closely to the approximate result in Appendix B. 

The particle energy in the transformed system is a constant, with the 
initial velocity given by the transformation velocity ( in the simple case 
where the initial velocities are approximately zero in the stationary system) . 
It is therefore apparent that protons and electrons will be ejected with the 
same velocity, a result which was found before (Appendix B) only for the 
first ana second approximations. 

An oscillation frequency about the neutral sheet is found approximately, 
(Equation (ap)) and the number of oscillations before ejection is determined 
Equation (aq) ) . It is seen that this number is proportional to the square 
root of the mass, so electrons will execute about l/40th the number of 
proton oscillations before ejection. 


- 7 - 



and 


Summary . Two models of possible current sheets are considered, 
approximate analytical results of charged particle trajectories about 
these sheets are found. 

In the simplest model, where the magnetic field varies linearly 
across a neutral sheet, and there is an electric field perpendicular to 
the magnetic field and parallel to the current sheet, charged particles of 
either sign execute damped oscillations about the sheet, accelerating 
along the sheet. Thus for this simple model, particles never come out of 
the sheet and their energies go to infinity. 

A new model is constructed by adding a small component of the 
magnetic field perpendicular to the sheet. The addition of this small 
component not only turns both protons and electrons toward the same 
direction 90 ° away from the accelerating electric field, but serves to 
eject the particles from the sheet. Protons and electrons are thus 
accelerated and are then ejected with the same velocity, the electrons 
being turned much faster and being ejected sooner than the protons. 

The pitch angle of ejected particles about a line of force is pro- 
portional to the size of the small perpendicular component cf magnetic 
field. Thus, if this component is small, all particles will be ejected 
nearly along lines of force . 

The energy gained is inversely proportional to the square of the 
perpendicular component of magnetic field. (See Equation (wj). Therefore 
in agreement with the simple model, the energies go to infinity when this 
component is zero. Moreover the energies of the ejected particles can be 
large if this component is small. 


- 7a - 



Appendix A . The first integrals of Equations (3) and (5) are: 


Z = 2 0 - C 3 1 



(a) 


and 


t (* 2 + 2 2 ) + C 3 z = A (x£ + Zq) + C 3 Z 0 

liquation (b) is just the equation of conservation of energy. The zero 
subscripted values refer to initial values. Equation ( 3 ) becomes, 
using (a): 


(b) 


x c, X Q-io + i C. (x^ - xp-r C 3 t) (c) 

For large enough time, the quantity in parenthesis on the right hand 
side of Equation (c) will be positive and monotonically increasing 
implying oscillatory, bounded motion in x(t). Thus for large time 
Equation (c) is approximately: 


The solution 


x « - C 1 C 3 x t 


q . - a b 

-V __ 

K m y d 


to Equation (d) is given by Jahnke and Erode [1945 > p p 147] 



( 

\ 


d) 


v e y 


- o 



and Z, is a linear combination of Bessel Functions of the First and 

73 

Second Kinds, of order one-third. Approximating (e) for large time, 
f jahnke and Em.de , 1945, PP • 133], 



A and B are constants depending on the initial values. 

z(t) may be obtained as a function of time by integrating Equation (a): 


z 


(t) 



(g) 


+ smaller oscillatory terms 


From (f), after large time the amplitude of oscillation decays as 
1 / 1 and from (g) and (b) it is seen that the kinetic energy increases 
as t . 

Appendix B . First integrals of Equations (8), (9), and (10) are 
easily 'obtained. They are: 


U 2 - i 2 ) = 


C 3 (z - z 0 ) - ~ U 2 - Zq) - c 2 Ti ! y 2 dt 


f , \ 


y - y 0 = c z n (z - z o) 


c. 


(z - Z n ) 


(x 2 - x 2 ) - C 3 Tl (y - y 0 ) 


- t 2 ^ A 0 ^ 


r . \ 

U/ 
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Equation (j) becorr.es, on another integration: 




'0 > ~ 


T J o x " dt 


c T| I ydt + C t 
2 J 0 4 


whe j 


C ,*o 


C 4 - zo + — j- + C 2 T]y 0 


(k) 


Using Equation (k) with Equation (i), we have: 


C -, ^ t- ^ 

- c 5 t 2 + C 6 t - c 7 I X 2 dt - c 2 -n ; ydt N - 


' 2 " J 0 ^ J 


( 1 ) 


wnere 


c 2 c 3 


C 1 C 2 


c 5 - 2 ’ ° 6 “ C 2 C 4 ’ c 7 - 2 


From Equation (9) we have: 


c 2 T| i y z dt = - (y 2 - y 2N ' 


(m) 


So that Equation (h) becomes the energy integral: 


y 2 + z 2 


z 2 ) + C 3 z = - (xq + + Zq) + C 3 z = 


E 0 


IB 


uO 


Using Equation (j), Equation (8) becomes: 


x = -k(t, x, y, Tl) 


(o) 


wnere 


r 2 


k(t, x, y, T|) = C t + x— 4- C 7)y + Cb 


(p) 


and 


C 8 “ ‘ 2 


C 2 y 0 - = - c 4 
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k goes to zero when 


t 


1 / C i x 2 

E7 + k ” y 



(q) 


r is therefore the critical time when k goes to zero and then becomes 
negative ejecting the particle. (z becomes positive.) 


Approxitr.ati ons . Integrating Equation (l ) , we obtain: 


C 6 t 2 Ti C 5 t 3 

y = y^ + yot + — 7 - - — 7 — ( r ) 

where the integral over x 2 in (l) is assumed small since x(t) is 

oscillating, and the integral over y in (i) is multiplied by Ti 2 , so will 
also be neglected in this first approximation. To facilitate finding 
the critical time of ejection t, initial conditions are chosen such that 
y 0 = y 0 = 0 ana zo = -C^x 2 /2, which implies that = C 6 = Cs = 0. Using 
Equations (q) and (r) with these initial conditions, t becomes: 


/6 



T|C 5 t 3 
y(t) = - — t — 

y(t) = - T:C 5 t 2 


(s) 


:o 


(u) 
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At t = T, Equations (t) and (u) become: 


y(- 


/6 a 



and 


ba 



Thus , in this first approximation, the ejection velocity is inde- 
pendent of (q/m) because electrons are ejected must sooner than protons. 
(As the next approximation, if y(t) from Equation (t) is used in the 
integral over y from Equation (k) it is found that T is increased by 
about 30ji , | y ( '•' ) 1 i s decreased by about 50$, and the ejection velocity 
is still independent of (q/m). The first approximation should therefore 
give the right order of magnitude.) 

Pitch angle distributions of emergent particles . In order to calcu- 
late the pitch angle of an emergent particle, we must estimate the maxi- 
mum value x can have at the time of ejection. Since x(t> oscillate- 

intii t = T, it is reasonable to assume that x « x n . Using this 

max u 

tssumption and Equation (w) with Equation \7) for the ir.ag.._tic field 
' ith x = d, we have : 


V = x-. e 

r~> J X 


Tib 


e 

y 


and the cosine of the pitch angle <y> becomes: 

■N- 

(B-V) 2*n^0 - 3 J 

cos a = —— — = 1 - u * 

dv — 

;o 


where u = a/’o, the bulk flow velocity exterior to the neutral sheet, and 
only terms of order Tj 2 are kept. For small a, cos a ~ j - a 2 12, so 


. i o _ 







j i 


(z) 


■ne nuziDsr s in nnuauion 


ould be decreased to 3/2 by the second appreciation mentioned above 


(z) cones from the 3 iu Equation (w) . 


This 


Aotendix C. 


norentz uransi creation can oe 


oiie i!_c luS 


given by Equations (2) and (7) to a system where the electric field is zero. 
This is possible for the nresent model because of the form of the fields 


and because B r (Equation (7)) is a constant. If the transformation velocity 
x 

is chosen as: 


v 


-a 

T)b 


A 


y 


(aa) 


she fields become in the transformed 


= 0 


h 


b A /, 


1 - v 2 /. 


.2 A 


X 

d 


y / 


( 8 - 0 / 


; the primes refer to the transformed system. 

The equations of motion corresponding to Equations ( 3 ) , ( 9 ) a 


LO) 


= C x'z 


(ac) 
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- c 2 T]'s 


( 8-d) 


-C, x'x' - Co7| 'y ! 


(ae) 


v!lere ^ ' = It -A - P/c 2 . 

If we borrow the previous result that the x co-ordinate oscillates 
until t — t > then the term -C^x'x' in Equation (ae) will be assumed small 
initially. Equations (ad) and (ae) then become: 


y' 




^ a: 


f) 


-wy’ 


[as) 


mere - caT) ‘ • rnese are jusr tne coupled equations for simple harmonic 
lotion, and the solutions are: 


y 


• ' . . - / 

yo cos cpc 1 - z Q sin urn • 


! = z 0 c ° 3 tut 1 


(ai) 


~ i 


cceroarec. no 


ror transiornmcion velocities 
the previous boundary conditions (Appendix B) be com 

Sjc> bcin 




m "G 110 1GX*Gs.IlSX ChT^lSG. 


yo = 


pb 3 


z 0 = 


— •[ Y, 


o, -'-o/2 


(aj) 


li- 



i* or 


anp; 


■I-’ y° can De J - ar o e compared to z Q , so Equations (ah) and (ai) 


y' 


r O COE 


f‘.-b f 


(oiO 


yo sin 1 


a) 


implying a circular drift in the 
drove negatively until cut' = rr/2, 

c o zero at Got * — tt. 


1 starts at zero and 
men decreases in absolute value 
positive thereafter, and so by the 


saroe arguments in the previous section and by Equation (ac) ve see tin 
mis nine is the ejection time. 


- v 


GO 


(q/m) bp 




°r velocities small compared to the speed of light, this uin 
he same m the two systems . me second approximation mentis 
replaced the /g in Equation ( s) by/ 




j_q vmcn is very 


Jl CiS jSO u. 


ty Ltci O .1. 


~ o^- v uo/ * v ren aquation (as) 


J St the time of e,iect_c 


- GU~ 


close uO . i 






tr 




^ ‘ J ~ ~ 7o* transforming back to the 


unprimed system gives 


y(T) = - 


vo 


(an) 


iac u or of / (1st approximation Appendix 3) or 3/2 (i 
pproximation Appendix E) . 



.-lore inf oration can. now be cotainea tnan an Appenuix B , oy 


anroxira'Cia’: z 


Equation (al) 


- yo /2 




Using Equation (ac) with Equation (ac) we see th 
oscillate with frequency 


co- ore 


will 


tu 


f c ^°\% fisMi) 1/2 

^ 2 j I 2yd / 


zp) 


nd the number of oscillations about the neutral sheet 


is eyeccea 




n 2b ^ 2(q7m)rj 3 dj 

^ J 


V/2 


lire above results can be extended to larger velocities (l; 
not mailing the approximation v « and keeping all of the 
; lorentz transformation. 
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